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' Recent computer simulation results [Barrat et ai, Physica A 334 (2004) 513] for granular mixtures 

. subject to stochastic driving have shown the validity of the Einstein relation e = D/{ToX) = 1 

' between the diffusion D and mobility A coefficients when the temperature of the gas T is replaced 

04 ' by the temperature of the impurity To in the usual Einstein relation. This problem is analyzed 

, in this paper by solving analytically the Boltzmann-Lorentz equation from the Chapman-Enskog 

CL|i method. The gas is heated by the action of an external driving force (thermostat) which does work 

' to compensate for the coUisional loss of energy. Two types of thermostats are considered: (a) a 

deterministic force proportional to the particle velocity (Gaussian thermostat), and (b) a white 
noise external force (stochastic thermostat). The diffusion and mobility coefficients are given in 
terms of the solutions of two linear integral equations, which are approximately solved up to the 
' (— I ' I second order in a Sonine polynomial expansion. The results show that the violation of the Einstein 

. relation (e ^ 1) is only due to the non-Maxwellian behavior of the impurity velocity distribution 

' function (absence of the Gibbs state). At a quantitative level, the kinetic theory results also show 

Ch ' that the deviation of e from 1 is more significant in the case of the Gaussian thermostat than in the 

case of the stochastic one, in which case the deviation of the Einstein relation is in general smaller 
than 1%. This conclusion agrees quite well with the results found in computer simulations. 
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I , I. INTRODUCTION 

Q . In spite of the growing theoretical support for the validity of the hydrodynamic description to granular fluids, 
^ ' some care is warranted in translating properties of normal fluids to those with inelastic collisions. Thus, for instance, 
in the case of an impurity particle immersed in an unforced granular gas the response to an external force on the 
impurity particle leads to a mobility coefficient that violates the Einstein relation between the mobility and diffusion 
^ , coefficients. This violation has been recently studied by Dufty and Garzo by using linear response methods and 

■ without any limitation on the mass ratio, density, or degree of inelasticity. The analysis shows that there are three 
y—i ' different reasons why the Einstein relation does not apply in granular gases: the absence of the Gibbs state, the cooling 

, of the reference state, and the occurrence of different kinetic temperatures between the impurity and gas particles. 

■ However, very recently Barrat et al. Q have performed computer simulations to test the validity of fluctuation- 
dissipation relations for a mixture constituted by two granular gases with different masses and restitution coefficients. 
The system is driven by an external energy input to achieve a stationary state. Specifically, they use a homogeneous 
stochastic driving (thermostat) in the form of a white noise plus a friction term. Given that the partial granular 

■ temperatures of each species are different, they observe that the classical equilibrium Green-Kubo relations are satisfied 
separately for each component, the role of the "equilibrium" temperature being played by the granular temperature 
of each component. In particular, in the limit in which the mole fraction of one of the species is negligible (tracer 

''O limit), the Einstein relation holds when the temperature of the gas T is replaced by the temperature of the impurity 
Ch Tq. In this case, the (modified) Einstein ratio between the diffusion coefficient D and the mobility coefficient A is 
verified, namely. 
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.- — = 1. (1) 

The results obtained by Barrat et al. ||3| from computer simulations motivate this work. The aim of this paper 
is to determine the Einstein ratio e for a driven granular gas in the context of the Boltzmann equation. As in the 
case of the free cooling gas studied in Ref. 0, I will use kinetic theory tools to explicitly get the dependence of e on 
the restitution coefficients as well as on the mechanical parameters of the system (masses and sizes). The theoretical 
estimates derived in this paper show that e 7^ 1, although the deviations of the Einstein ratio from unity are in general 
quite small. This conclusion is in good agreement with the observations made by Barrat et al. 
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The paper is organized as follows. Section [H] describes the problem I am interested in. In particular, two types of 
external forces (thermostats) to inject energy to the gas and reach a steady state are considered: the Gaussian and 
the stochastic thermostats. In the case of the Gaussian thermostat, the gas is heated by the action of an external force 
proportional to the velocity. This type of "anti-drag" force can be justified by Gauss's principle of least constraint and 
has been widely used in noncquilibrium molecular dynamics simulations of normal fluids 4, 5] . Another mechanism of 
heating (similar to the one used in Ref. Q) is to consider the action of a random external force, which gives frequent 
kicks to each particle between collisions. If this stochastic force has the properties of a white noise, it gives rise to a 
Fokker-Planck diffusion term in the Boltzmann equation Section llTTl is devoted to the calculation of the diffusion D 
and mobility A coefficients by solving the Boltzmann-Lorentz equation by means of the Chapman-Enskog method . 
Both coefficients are explicitly determined by using both types of thermostats, being given in terms of the solutions 
of two linear integral equations. A practical evaluation of both transport coefficients is possible by using a Sonine 
polynomial expansion. Here, D and A have been computed in the second Sonine approximation (namely, retaining 
two Sonine polynomials). As said above, in the case of the stochastic thermostat, the analysis agrees with the results 
of Barrat et al. Q- In the case of the Gaussian thermostat, the deviations of e from unity are more significant, 
although they are still small for not too inelastic systems. The paper is closed in Section Hvl with a brief discussion 
on the results obtained in this paper. 



II. DESCRIPTION OF THE PROBLEM 



Let us consider a granular gas of smooth inelastic hard spheres (of mass m, diameter cr, and interparticle coefficient 
of restitution a) in a homogeneous state. In the low-density regime, its velocity distribution function /(v) obeys the 
nonlinear Boltzmann equation Due to dissipation in collisions, the gas is cooling unless a mechanism of energy 
input is externally introduced to compensate for coUisional cooling. In experiments the energy is typically injected 
through the boundaries yielding a inhomogeneous steady state. To avoid the complication of strong temperature 
heterogeneities, it is usual to consider the action of homogeneous external (driving) forces acting locally on each 
particle. These forces are called thermostats and depend on the state of the system. In this situation, the steady-state 
Boltzmann equation reads 

^/(v)- J[v|/,/], (2) 
where the Boltzmann collision operator J[v|/, /] is 

J [Vll/, /] - / j da e(^ • gi2)(^ • gl2) [«-'/(v'i)/(v^) - /(Vi)/(V2)] . (3) 

In Eq. is an operator representing the effect of the external force, ? is a unit vector along the line of centers 

of the two colliding particles, Q is the Heaviside step function, and gi2 = Vi — V2. In addition, the primes on the 
velocities denote the initial values {v5^,V2} that lead to {vi,V2} following a binary collision: 

v'l = vi - i (l + (S • gi2)5^, V2 = V2 + ^ (l + a"^) (o' • gi2)5^ . (4) 

Several types of thermostats can be used. As said in the Introduction, here I will consider two simple choices. One 
of them is a deterministic thermostat based on Gauss's principle of least constraints ^|^. In this case, T has the 
form 

•^/(v) = ■ [v/(v)] , (5) 



where C, is the cooling rate due to collisions 



2,nT 

Here, m is the mass of a particle and T is the granular temperature defined as 



C = -^ i dvmv'JWJ]. (6) 



n being the number density 



\nT= I rfv|.V(v), (7) 



dv/(v). (8) 
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FIG. 1: Plot of the coefficient co versus tfie coefBcient of restitution a 
thermostat (a) and the stochastic thermostat (b). 



Qo for mo/m — 0.5 and ao/a — 1 for the Gaussian 



Another way of heating the gas is by the action of a random force in the form of an uncorrelated white noise 9] . The 
corresponding operator has a Fokker-Planck form 



2 m V av 



/(v). 



(9) 



The use of this stochastic force has attracted the attention of many theorists in the past years to study different prob- 
lems, such as non-Gaussian properties (cumulants, high-energy tails) IE [13, 

long-range correlations lUj, coUisional 
statistics and short-scale structure ^3 : transport properties [13 , and fluctuation-dissipation relations jl4l| . 

So far, the exact solution to the Boltzmann equation |(2J) is not known, although a good approximation for thermal 
velocities can be obtained from an expansion in Sonine polynomials. In the leading order, / is given by 
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(10) 



where v* = v/vth, Vth = y^iT/m being the thermal velocity. The coefficient c, which measures the deviation of / 
from the Maxwellian form, is related to the kurtosis of the distribution. Its value depends on the thermostat used 
and has been estimated from the Boltzmann equation up to linear order in c. In the case of the Gaussian thermostat 
the coefficient c is given by 



c(a) 



32(l-a)(l-2a2) 
81-17a + 30a2(l-a)' 



(11) 



Estimate (jllf) presents quite a good agreement with Monte Carlo simulations of the Boltzmann equation [Tol ITsIl . 
It is interesting to remark that in the homogenous problem the results obtained with and without the Gaussian 
thermostat are completely equivalent when one scales the particle velocity with respect to the thermal velocity [T(tI |. 
As a consequence, result (|llll applies to the free cooling case as well. In the case of the stochastic thermostat © the 
result is 



c{a) 



32{1 - a)il - 2a^) 



(12) 



241 - 177a + 30^2(2 -a)' 

As in the Gaussian case, Eq. (|12() agrees quite well with computer simulations In general, both theoretical and 
numerical calculations show that, for not too inelastic systems, c is very small. Deviations from the Maxwellian 
form are important, however, if one considers the tails of the distribution functioU j w hich are exponential for the 
stochastic thermostat and a "stretched" Gaussian for the Gaussian thermostat [a, [la ■ The cooling rate ^ can also 
be determined from the Sonine approximation H10|l with the result 



1 + A, 
32 



(13) 
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FIG. 2: Plot of the coefficient co versus tfie coefBcient of restitution a = ao for mo/m — 5 and ctq/o" = 1 for the Gaussian 
thermostat (a) and the stochastic thermostat (b). 



We assume now that a few impurities (of mass mg and diameter ctq) are added to the system. Given that their 
molar fraction is negligible, the state of the granular gas is not disturbed by the presence of impurities and so the 
velocity distribution function /(v) obeys the Boltzmann equation Moreover, one can also neglect collisions among 
impurities themselves versus the impurity-gas collisions, which are characterized by a coefficient of restitution ao- We 
want to analyze the diffusion of impurities immersed in the granular gas when the current of impurities is generated 
by the presence of a weak concentration gradient Vng and/or a weak external field E (e.g. gravity or an electric 
field) acting only on the impurity particles. Under these conditions, the velocity distribution function /o(r,v, i) of 
impurities verifies the Boltzmann-Lorcntz equation 

dtfo + V • V/o + — • ^/o + ^/o = J[v|/o, /], (14) 
mo av 



where 



JWllfoJ] = j d^2 j rfSe(CT • gi2)(CT • gi2) 

X [aoVo(r,v'i;t)/(r,v^;i)-/o(r,vi;t)/(r,V2;0] , (15) 



a = {a + o'o)/2, and 



Vi=vi ■ (1 + ao )('^-gi2)cr, V2=V2H ■ (1 + ) cr • gi2)cr . (16) 

m + mo ^ " ' rn + niQ ^ ' 

The partial temperature of impurities Tq is defined as 

^noTo = 1 dvI|I^;Vo(v), (17) 
where no is the number density of impurities, 

«o = y dv/o(v). (18) 
In addition, it is convenient to introduce the cooling rate Co associated with the partial temperature To of impurities: 



Co = -1^^ \ dvmov'J[v\foJ]. (19) 
3noTo J 

Momentum and energy are not collisional invariants of the Boltzmann-Lorentz collision operator J[fo, /]. Only the 
number density of impurities is conserved: 

dtno + \/-h = 0, (20) 
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FIG. 3: Plot of the coefficient co versus the mass ratio mo/m for a = ao = 0.8 and ao/a = 1 for the Gaussian thermostat (a) 
and the stochastic thermostat (b). 



where the current of impurities jo is defined as 



Jo 



dvv/o(v). 



(21) 



The conservation equation H20() becomes a closed hydrodynamic equation for no once jo is expressed as a functional 
of no and the external field E. 

In the absence of diffusion (homogeneous steady state), Eq. H14|) becomes 



^/o = J[v|/o,/]. 



(22) 



This equation has been recently analyzed by using both types of thermostats 0, 0| . The results show that the 
temperatures of the gas and impurities are clearly different (T ^ Tq) and so, the energy equipartition is broken down. 
The failure of the energy equipartition in granular gases has been confirmed by computer simulations |l8l l20l I2H 
and even observed in real experiments of vibrated mixtures in two and three |23| dimensions. In general, the 
temperature ratio 7 = Tq/T presents a complex dependence on the parameters of the problem. However, according 
to Eqs. (jSJ,©, and (|22|) the condition for determining the temperature ratio is different for each type of thermostat. 
In the case of the Gaussian thermostat (jSJ, when one multiplies both sides of Eq. (|22|l by and integrates over v, 
the steady state condition yields the equality of the cooling rates 0, |3 



In the case of the stochastic thermostat 0, 7 is obtained from the condition 

moC = TO7Co- 



(23) 



(24) 



Requirements (|23|) and 124(1 lead to a different dependence of the temperature ratio Tq/T on the control parameters, 
namely, the mass ratio mo/rn , the size ratio cto/ct and the coefficients of restitution a and ao- 
A good estimate of 7 can be obtained by considering the first Sonine approximation to /q: 



/o(v) -> non'^'v-^'ey^e 



Co 



15 

T 



(25) 



where = moT/mTo is the mean square velocity of the gas particles relative to that of impurities. The coefficient co 
can be determined by substitution of Eqs. 1)10(1 and ((25(1 into the Boltzmann-Lorentz equation 1(22(1 and retaining all 
terms linear in cq and c. The explicit expression of the coefficient cq is given in Appendix ^ Once the coefficient co 
is known, the temperature ratio can be obtained from the constraint ((23(1 for the Gaussian thermostat and II24I) for 
the stochastic thermostat. 

Before illustrating the dependence of Cq and 7 on the coefficients of restitution and the parameters of the mixture, 
it is instructive to consider some special limit cases. First, when impurities are mechanically equivalent to the gas 
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FIG. 4: Plot of the temperature ratio Tq/T versus the coefficient of restitution a = ao for mo/m — 0.5 and ctq/o" ~ f for the 
Gaussian thermostat (a) and the stochastic thermostat (b). 



particles (m = mo, a = ao, ckq = a), the results of the single gas case Q are recovered, namely, T — To, and 



32(l-ao)(l-2ag) 
81 - 17ao + 30a^(l - ao)' 



for the Gaussian thermostat while 



^ 32(l-ao)(l-2ag) 

° 241 - 177ao + 30a^(l-ao)' ^ ' 

for the stochastic thermostat. Also, when the collisions among gas particles are elastic (a = 1 and so c = 0), our 
Sonine solution gives Cq = and 

^= 1±^2 (28) 

2 + (1 - Q;o)(TO/mo) 

for both types of thermostats. This result agrees with the exact solution of the Boltzmann-Lorentz equation found 
by Martin and Piasecki in the study of the stationary state of a test particle immersed in a homogeneous fluid in 
equilibrium, undergoing dissipative collisions with the fluid particles. Finally, in the weak dissipation limit (1 — a ^ 1, 
1 — ao ^ 1)) the temperature ratio 7 behaves as 

^^l-«), (29) 
for the Gaussian thermostat and 




m + mo, . (ay i /mo(m + moj\ , , , 

for the stochastic thermostat. The two latter terms on the right hand side of Eqs. H29|) and H3U|I represent two different 
types of inelastic collisions providing independent mechanisms to enforce the breakdown of energy equipartition. 

Beyond the above limit cases, it is quite difficult to provide the explicit dependence of cq and 7 on the characteristic 
parameters of the system. For the sake of concreteness, now we consider the case a = ao- Figures Q and |21 show the 
coefficient co as a function of ao for ao/a — 1 and the mass ratios mo/m — 0.5 and mo/m = 5, respectively. We 
consider the cases of the Gaussian and stochastic thermostats. When the mass ratio is smaller than 1, we see that 
the dependence of cq on ao is qualitatively similar in both thermostats, its value being small for not too inelastic 
systems. This means that for moderate dissipation the distribution function fo of the homogeneous state is quite 
close to a Maxwellian at the temperature of the impurity particle To. However, according to Fig.|21 the dependence of 
Co on ao is quantitatively different for both thermostats when impurities are heavier than the particles of the gas. In 
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FIG. 5: Plot of the temperature ratio Tq/T versus the coefficient of restitution a = qq for mo/m = 5 and ao/a — 1 for the 
Gaussian thermostat (a) and the stochastic thermostat (b). 

particular, the curves of Fig. [21 exhibit a different convexity in the quasielastic limit. To know something more on the 
dependence of cq on the mass ratio, in Fig.|51we plot cq versus the mass ratio nio/rn for a = ao = 0.8 and cro/a = 1. 
The qualitative dependence of cq on mo/m is quite similar for both thermostats since first cq increases with mo/m, 
reaches a maximum and then decreases with the mass ratio. At a quantitative level, we observe that the magnitude 
of Co for the stochastic thermostat is in general smaller than that of the Gaussian thermostat and consequently, the 
deviation of /o from its Maxwellian form is more significant in the case of the Gaussian thermostat than in the case 
of the stochastic force. The dependence of the temperature ratio Tq/T on the restitution coefficient ao is illustrated 
in Figs. ^ and |S1 for the cases mo/m = 0.5 and cto/ct = 1 and rrio/m = 5 and (Jq/cf = 1, respectively. We consider 
again both types of thermostats. We observe that the temperature of the impurity is smaller (larger) than that of 
the gas when the impurity is lighter (heavier) than the particles of the gas, whatever the thermostat considered is. 
The lack of energy equipartition is quite apparent in all the cases studied, especially when the mass ratio mo/m is 
larger than 1 in the case of the Gaussian thermostat. This is indicative of what happens in the homogeneous cooling 
state in the large impurity/gas mass ratio where a peculiar "phase transition" has been recently found for which the 
diffusion coefficient is normal in one phase but grows without bound in the other |26|. It must be remarked that the 
predictions obtained for cq and 7 from the leading Sonine approach 125|) compare quite well with recent computer 
simulations |19L.2Qj. 

III. DIFFUSION AND MOBILITY COEFFICIENTS 

As said in the Introduction, my aim is to get the diffusion and mobility coefficients in the limit of small concentration 
gradient Vng and weak external fields. These transport coefficients are determined by solving the Boltzmann-Lorentz 
equation (|14|l from the Chapman-Enskog method 0|. This method assumes the existence of a normal solution in 
which all the space and time dependence of /o is through the hydrodynamic fields. In this problem, the normal 
solution is explicitly generated by expanding /o in powers of the gradient VriQ and the field strength E; 

/o-/r+^/.i'^+--- , (31) 

where each factor S corresponds to the implicit factors VriQ and E. The time derivative is also expanded as dt ~ 
9f + S d^^^ + • ■ • , where the action of the different operators d^.''^ can be obtained from the balance equation H20|l . 
They are given by 

di"'>no^O, d^,'^no = -V-j^t'\ k>l, (32) 

where 




(33) 
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The zeroth-order approximation f^^^ is the sohition of Eq. (|22|l , whose approximate form is given by Eq. ^ but 



taking into account now the local dependence on the density uq. Since /q"'' is isotropic, it follows that the flux of 

ler, i.e., Jq^'' = 0, anc 



impurities vanishes at this order, i.e., jl^^ = 0, and so d^^^^no — 0. To first order in 6, one has the kinetic equation 



E d 



(0) 



V . V + — • — 

mo ov / 

^ E d' 

v(Vln7io) + -TT- 

niQ ov 



^0 



(0) 



(34) 



The second equality follows from the space dependence of flf^ through no- The solution to Eq. (|34|) has the form 



/^^^(v) =^(v). Vlnno + B(v).E. 
Substitution of Eq. H35(l into Eq. (|34|l leads to the set of linear integral equations 

TA-J[v\AJ] = -vfl,°\ 



(35) 



(36) 



Too ov 

To first order, the current has the structure 

ii^'' = -L>Vlnno + AE, 



(37) 



(38) 



where D is the diffusion coefficient and A is the mobility coefficient. According to Eqs. (|33|l and H35|) . these coefficients 
are given by 

D = -^[d^^-A, (39) 



A = - i dwM -B. 

In the elastic limit (a = ao = 1), T = Toj a-nd so /q"'' is the local equilibrium distribution, i.e 



.(0). ^ / "lo / 



mQV 
2T 



(40) 



(41) 



In this case, dfl^^ /dw = —{mi:,w/T)fl^'^ and the integral equations l|36(l and (|37|l yield A - 
one recovers the usual Einstein relation between the diffusion and mobility coefficients 

D 



TB. As a consequence. 



(42) 



However, at finite inelasticity the relationship between D and A is no longer simple. In the driven case, there 
are two separate reasons for which the Einstein relation in its usual form is violated. First, since /q''^ is not a 
Gaussian distribution then, dfl^^ /d\' ^ —{mQv/T)/^^'' and D is not proportional to A. The second reason is that 
the temperature of the impurity particle is different from that of the fluid, except when the impurity is mechanically 
equivalent to particles of the fluid. This last source of discrepancy could be perhaps eliminated if one replaces the 
temperature of the gas T by the temperature of impurities To in the usual Einstein ratio (|42|l . This change leads 
to the modified Einstein ratio given by Eq. QJ. As a matter of fact, if one took the Maxwellian approximation l|41(l 

for flf ^ (with To instead of T), then dfj^^ /dw — — (toov/To)/o*''' and the modified Einstein relation between the 
diffusion and mobility coefficients would apply (e = 1). Consequently, the only reason for which e ^ 1 is due to 

the non-Maxwellian behavior of the velocity distribution function (absence of the Gibbs state), which leads to 
small deviations that could be difficult to detect in computer simulation experiments. This conclusion agrees with 
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the results obtained from recent molecular dynamics simulations [3, where no deviations from the Einstein relation 
e = 1 were observed for a wide range of values of the coefficients of restitution and parameters of the system. 

For practical purposes, the linear integral equations and 1)37(1 can be solved by using the first terms in a Sonine 
polynomial expansion. The reliability of this Sonine approximation has been confirmed in several inhomogeneous 
situations. Thus, in the case of the shear viscosity coefhcient for a single gas, the first Sonine predictions JJi, 2^ 
compare quite well with a numerical solution of the Boltzmann equation in the uniform shear flow state obtained 
from the direct simulation Monte Carlo (DSMC) method Such a good agreement for the shear viscosity has been 
also extended to the case of multicomponent systems in the low-density limit and for finite densities S^j . More 
recently, the tracer diffusion coefficient has been evaluated up to the second order of the Sonine expansion [13 . As in 
the case of the shear viscosity coefficient, the Sonine solution agrees quite well with computer simulation results over 
a wide region of the parameter space. In this paper, the diffusion and mobility coefficients will be determined up to 
the second Sonine approximation. In this case, the unknowns ^ and 3 are approximated by 

A{v) ^ /o,m(v) [aiv + a2So(v)] , (43) 



e(v) ^ /o,m(v) [61 V + 62So(v)] , (44) 
where /o,m(v) is a Maxwellian distribution at the temperature Tq of the impurities, i.e. 



and So(v) is the polynomial 



/o,M(v)=no(^) cxp(-^), (45) 



2 " 2 



The coefficients ai, 61, 02, and 62 are defined as 



™« ^rfvv.r^l, (47) 



hi J SnoTo J \B 



' ^d.So-(i]. (48) 



b2 J 15noT3 J " \ B 

These coefficients are determined by substitution of Eqs. H43|) and (|44|l into the integral equations (|36() and 137|l . 
respectively. The details are carried out in Appendix In the case of the Gaussian thermostat, the Einstein ratio 
(P) becomes 

(49) 

where C*, 1^2 and 1/4 are explicitly given in Appendix IbI In the case of the stochastic thermostat, one gets 

. = 1-1^. (50) 



Equations H49I) and (|50|l clearly show that the deviations from the Einstein relation e = 1 are proportional to the 
coefficient cg, which measures the deviation of the velocity distribution function of impurities from the Maxwellian 
distribution. Both c and cq vanish in the limit of elastic collisions and otherwise give contributions of the order of a 
few percent to 1^2 j and C*, even at very strong dissipation. According to Eqs. (|49|l and (|50|l it is apparent that the 
Einstein ratio exhibits a complex dependence on a, ao, and the mechanical parameters of the mixture (the mass and 
size ratios). A simpler form for the Einstein ratio corresponds to the case of mechanically equivalent particles (i.e., 
m = mo, cr = CToi and a = ao). In this case. To = T, cq = c, and Eqs. (|49|l and H5Ufl reduce in the quasielastic limit 
to a common ratio given by 

^-1 + TT^(1-")- (51) 
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TABLE L Values of the Einstein ratio e for the different cases studied in Ref. Q. The last two columns refer to the results 
obtained by using the Gaussian and stochastic thermostats, respectively. 



e 




«0 

FIG. 6: Plot of the Einstein ratio e versus the coefficient of restitution qq for the Gaussian thermostat in the cases: (a) q = ao, 
mo/m = 5 and cto/it — 1; (b) a = qq, mo/m — 0.5 and ao/a = 1; and (c) a = 0.5, mo/m = 10 and ao/a — 1. 

Equation l|51|l shows most directly the dependence of e on dissipation. Also, when the gas is at equilibrium (a — 1), 
then cq — with 7 given by Eq. (|28|l and so, the Einstein relation is exactly verified. 

As said in the Introduction, very recently Barrat et al. have used the DSMC method [2^ and molecular dynamics 
simulations to investigate the validity of fluctuation-dissipation relations for a binary mixture of inelastic hard disks 
with identical diameters but different physical properties. The system is also subjected to a external force which has 
the form of a white noise [Eq. ©] with the addition of a viscous term. In the single tracer case, their computer 
simulation results clearly support the Einstein relation Q. Although a quantitative comparison between the theory 
presented here and their results is not possible, it is illustrative to get the values of the Einstein ratio e predicted by 
the theory for the same cases as considered in the simulations. The results are displayed in Tabled for the Gaussian 
and stochastic thermostats. Note that the latter coincides with the thermostat used in simulations, except for the 
viscous term. We see that the deviations from the Einstein relation are very small (less than 1%), even for moderate 
dissipation (say for instance, a ~ ao — 0.7). This confirms the conclusions obtained in Ref. Q from computer 
simulations. 

To illustrate the influence of dissipation on the Einstein ratio more generally, in Figs.Eland[7|the Einstein ratio e is 
plotted versus ao for ao/a — 1 and different values of the mass ratio mo/m. Figure El shows the results obtained by 
using the Gaussian thermostat while Fig. [7| refers to the results obtained from the stochastic thermostat. Although 
the mechanical properties between both species (impurities and particles of the gas) are quite different, we observe 
that e is very close to 1, especially in the case of the stochastic thermostat where again the deviations are smaller 
than 1%. However, in the case of the Gaussian thermostat the deviations from unity are even about 8% (for instance, 
for a ^ ao = 0.5 and mo/m = 10). This failure of the Einstein relation could be detected in computer simulations by 
using the Gaussian thermostat instead of stochastic driving. In addition. Figs. and [7| illustrate again the fact that 
the transport properties are affected by the thermostat introduced so that the latter does not play a neutral role in 
the problem [l3l| . 
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FIG. 7: Plot of the Einstein ratio e versus the coefficient of restitution ao for the stochastic thermostat in the cases: (a) a — uq, 
mo/m — 5 and ao/cr = 1; (b) a — ao, mo/m = 0.5 and ao/a — 1; and (c) a — 0.5, mo/m — 10 and ao/cr = 1. 

IV. DISCUSSION 

The goal of this paper has been to analyze the validity of the Einstein relation between the diffusion and mobility 
coefficients in a heated granular gas. This analysis has been mainly motivated by recent computer simulations 
performed by Barrat et al. ^] where no apparent deviations of the Einstein ratio from unity were found when the 
temperature of the gas T is replaced by the temperature of impurities Tq in the usual Einstein ratio, Eq. This 
conclusion contrasts with the one obtained before by Dufty and Garzo in the unforced gas case where it was clearly 
shown the violation of the Einstein relation. The analysis carried out in Ref. indicates that the deviations of the 
Einstein ratio in its usual form from unity has three distinct origins: non-Gaussianity of the distribution function of 
the homogeneous cooling state, energy nonequipartition, and time evolution of the granular temperature. Given that 
the last source of discrepancy is no present in the driven case, it is interesting to assess the effect on the violation of 
the Einstein relation coming from the first two reasons. 

To achieve a steady state, the granular gas is driven by the action of external forces (thermostats) that accelerate the 
particles and hence compensate for collisional cooling. Here, I have used two types of thermostats: (a) an "anti-drag" 
force proportional to the particle velocity (Gaussian force), and (b) a stochastic force, which gives frequent kicks to 
each particle between collisions. The introduction of these thermostats has the advantage of avoiding the intrinsic 
time dependence of the temperature in the free cooling state, but at the price of introducing new effects induced by the 
external forcing. The diffusion D and mobility A coefficients have been determined by solving the Boltzmann-Lorentz 
equation by means of the Chapman-Enskog method. In the first order of expansion, D and A are given in terms of 
the solutions of two linear integral equations. An approximate solution of these integral equations has been obtained 
by expanding the distribution function of impurities in a Sonine polynomial expansion. Here, I have retained two 
Sonine polynomials (second Sonine approximation). As happens for the transport properties of a single gas the 
Sonine solution shows that the dependence of D and A on the parameters of the problem is different for each type 
of thermostat. Although not widely recognized, this conclusion illustrates the fact that generally the inclusion of an 
external force depending on the state of the system changes the apparent transport coefficients. 

As Barrat et al. P| suggests, the influence of the temperature differences between impurities and particles of the gas 
on the Einstein ratio can be mitigated in part if one uses Tq instead of T in the usual Einstein ratio. For this reason, 
I have explored the validity of the modified Einstein ratio (Q). The analytical results (based on the second Sonine 
approximation) show that the deviation of the Einstein ratio from unity is only due to the factor cq which measures 
the deviations of the velocity distribution function of the impurities from a Maxwellian. In fact, if one neglects the 
small effect of this factor, then e ~ 1. In addition, the results provided here also show that the deviation of e from 1 is 
more significant in the case of the Gaussian thermostat than in the case of the stochastic one. This is due essentially 
to the fact that in general the magnitude of the coefficient cq of the stochastic thermostat is smaller than that of the 
Gaussian thermostat. As an illustration. Figs. Eland [7| show that the deviations of the Einstein ratio from unity are 
generally smaller than 1% for the stochastic thermostat while these deviations are even about 8% in the case of the 
Gaussian thermostat. Given that the simulations carried in Ref. Q were made by using a stochastic driving (with 
the addition of a viscous term), the results of this paper agree with their observations since deviations smaller than 
1% are quite difficult to detect in computer simulations. 



12 



The results obtained here have been restricted to the single tracer case, i.e., a especial limit case of a binary mixture. 
Barrat et aZ.'s computer simulations also show the validity of the Einstein relations separately for each species, each 
one with its own temperature. I plan to analyze this issue in the next future. 
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APPENDIX A: EXPLICIT EXPRESSIONS OF Co AND co 

In this Appendix, the expressions of the partial cooling rate Co and the coefficient cq in the homogenous states 
driven by the Gaussian and stochastic thermostats are given. Most of these results have been recently reported by 
the author in Ref. the homogeneous cooling state. By using the leading Sonine approximations (|10|l and H25|l 

and neglecting nonlinear terms in c and cq, Co can be written as [l9ll2^ 



Co = Aoo + Aoico + A02C, (Al) 



where 



3 



1 



1/2 



^00 = ^V^nW'vthfi ( ] (1 + ao) 1-^(1 + ao)(l + ^) , (A2) 



2 



Aoi = ^V^na^va,fi^-^-^^j^{l + «o) [2(3 + 49) - + «o)(l + d)] , (A3) 



\o2 = -j^V^na^vthfii—^j (l + ao)[2 + 3/i(l + ao)(l + e)]. (A4) 



Here, a — {a + o'o)/2, /i = m/(m + ttiq), and 9 = moT/mTo. 
The coefficient cq is defined by 



15 \4noT§ J •'° 4 



(A5) 



This coefficient is determined from the Boltzmann-Lorentz equation H22II by multiplying that equation by v'^, and 
integrating over the velocity. When only linear terms in co and c are retained, the result in the case of the Gaussian 
thermostat is found to be [23 

-y0Co(l + y) =f^oo + f^oiCo + r!o2C, (A6) 

while in the case of the stochastic thermostat the result is 

15 

-yPCo = f^oo + f^oiCo + f^o2C, (A7) 

where use has been made of the conditions H23|) and (|24|) . respectively. Here, ^0 ^ 1 — ^ = mo/{m + mo) and I have 
introduced the quantities [19| 

Ooo = 2V^Tla\th^J-l^l^^-^^^^{l + ao)[-2{6 + 59) + ^l{l + ao)il + 9){U + 59) 

-8a*' (1 + aof (1 + 9f + 2/i3 (1 + aof (1 + 9f] , (A8) 
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^01 = ^ng^^thA^gA^ ^^^^f/, (1 + ao) [-2 (90 + 2310 + 184^2 ^ 40^3^ 

+3m (1 + ao) (1 + 9) (70 + 1176' + 446*2) - 24^^ _^ ^^)2 _^ ^,^2 ^ 

+30^^^ (1 + ao)^ (1 + 61)^] , (A9) 

^02 = ^na^w^M ^^^gf/^ ^ (l + ao)[2(2 + 5g) + 3Af(l + ao)(l + g)(2 + 50) 

-24^2 (l + ao)^(l + 6l)V30/i^(l + ao)^ (1 + 61)^] . (AlO) 

The final expression of cq is obtained by substitution of Eq. (jAlj) into Eq. 1IA6P for the Gaussian thermostat and Eq. 
l|A7p for the stochastic thermostat and neglecting nonlinear terms in cq and c. The result is 

Aoo + A02C + ^A^o ^^'^ (^^00 + f^02c) 
= lA + 2 -2.20 ' 

2 Aqo + Aoi + i5A*0 ^^01 

for the Gaussian thermostat while 



Aoo + A 02C+ ^^/6i2 (»oo + r2o2c) 
Aoi + j^fJ-o ^^2^01 



for the stochastic thermostat. 

Once the coefficient cq is given in terms of 7 = fio/{fi9) and the parameters of the mixture, the temperature ratio 
7 can be explicitly obtained by numerically solving the condition (|23|l for the Gaussian thermostat or the condition 
(|24(l for the stochastic thermostat. 

APPENDIX B: SECOND SONINE APPROXIMATION 

In this Appendix the coefficients ai and bi in the second Sonine approximation are determined. The knowledge of 
these coefficients allows one to get the diffusion D and mobility A coefficients from Eqs. (|39ll . (|40ll . and (|47|l . Let us 
consider first the case of the Gaussian thermostat jSjl. For this thermostat, Eqs. (|36|l and (|37|) become, respectively 

ic|^-(v^)- j[A/]=-v/r, (Bi) 

• (vB) - J[BJ] = - — (B2) 
2 av mo ov 

Substitution of Eqs. (|43|l and 1441) into the above integral equations gives 

1 

2^av ' ^ ^^'^f^-M^ + a2/o,J\/So) - ai^[/o,AfV, /] - a2J[/o,AfSo, /] = -v/($°\ (B3) 

•v(6i/o.MV + 62/0,^80) - 6iJ[/o,A./v,/] - fe2J[/o,A./So,/] = - — (B4) 
Next, Eqs. (|lj3() and (|lj4|) are multiplied by v and integrated over the velocity. The result is 

iyi - + ^^202 = -1, (B5) 

{^1 - \Obi + y2b2 = T-\ (B6) 
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where I have introduced the quantities 

7/1 = — 

3noTo 



VI ^- I dv V • J[/o, A/ V, /] , (B7) 



1^2 - f dvv J[/o,A/So, /]. (B8) 



If only the first Sonine correction is retained (which means 02 = 62 = 0), the solution to Eqs. ljB5|l and ljB6p yields 
ai = —Tobi, and so the modified Einstein relation applies (e — s- 1). 

To close the problem in the second Sonine approximation, one multiplies Eqs. (|B3p and IIB4(I by So(v) and integrates 
over the velocity. Following identical mathematical steps as those made before, one gets 

- CTq-') ai + [>^4 - ^C) "2 = -lcoT,-\ (B9) 
- (To') h + (1^4 - 62 = 0, (BIO) 



where I have introduced the quantities 



;/4 



4^ / dvSo- J[/o,mSo,/]. (B12) 
15 nolo J 



In reduced units and by using matrix notation, Eqs. (|B5I) and l|B9(l . along with Eqs. ljB6|l and ljB10|l . can be 
rewritten, respectively, as 

(B13) 



(B14) 





2*5 




Here, C* = C/z/, i/J" = t^2 = i^2/Tqv, v'^ = Tqv^/v, and 1/4 = i^4/i^ with v ~ na^vth- Further, — vai, 

0-2 — TQi'a2, bl = T^vhi, and h*2 = Tovh2. Equation l|B13p coincides with the one recently derived in the free cooling 
case ,231 • The solution to Eqs. (|B13p and ljB14p provides the explicit expression of the Einstein ratio in the second 
Sonine approximation. The result is 

e=^ = -^ = l-^ . (B15) 

The analysis in the case of the stochastic thermostat Q is similar to that made above for the Gaussian one. The 
linear integral equations are now given by 

A-J[AJ] = -.f^^\ (B16) 

2m \ov J 
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The corresponding matrix equations defining the Sonine coefficients {01,02,61,62} can be obtained by multiplying 
Eqs. (|ljl6() and Hijl7|l by v and So(v). After some algebra, one gets 



"1 - 



Co* < 



Co/2 



(B18) 



V - Co* 4 j V V V ' 

where = Qq/v and use has been made of the relation 

Cr _ CqTo 
TO mo 

The solution to Eqs. (|B18(I and (|B19() leads to the following expression for the Einstein ratio: 

^ 1 _ 



(B19) 



(B20) 



(B21) 

The collision integrals i/| and appearing in the Einstein ratios (|B15|) and (|B21|) have been recently evaluated by 
the author 0. They are given by 



y*2 = lv^y^\-] (i + «o) [0(1 + 0)]''/' 



16 V 1 



(B22) 



where 



A 



(2af, ~ 3ao + 4) 



3/2 



A-5- 



-B 



-/i(l + 9) [2/361-2(8 + 56')(7ao - U) + 2e'"i(29ao - 37) - 25(1 - ao)] 
+18P^e~^{8 + 56) + 2/361-1(25 + 666*) + 59-\6 + IW) - 5(1 + 9)0-^(6 + 56 

+ ^ (1 + 9)-^ {150^(1 + «o)(4/i(l + ao) - 5) 
Id 

+2 [45 + 540^^0 + 16M(ao - 36) + 4/1^(134 + 5ao + 6ag) 
-4^o(148 + Ai(7ao - 263))] + 9^ [-30 - Ai(267 + 217ao) 
+14^2(^7 _^ 29ao + 120^) + 10/io(7/i(l + ao) - 5)] 
+9 [-315 + 270Aio - 2^(55ao + 57) + ^^(449 + 326ao + 156a^) 
+2/io(-2 + ^(7ao + 277))]}, 

B = 5(1 + 2/3) + /i(l + 9) [5(1 - ao) - 2(7ao - 11)P9-'^] + 18(3'^9-^ 
+2^2(2^2 _ _^ 4)61-1(1 + 61)2 - 561-1(1 + 9) 

^YeJlTW ^^^^^^^ + t^''^^ + ao) - 5] + [2^0 {7^^{l + ao) - 5) 
+^l (-5(9 + 7ao) + m(38 + 62ao + 24a2))] - 15 + bAfil - 20/i(3 + ao) 
+2^2(40 ^ i9q,„ + Qal) + 2^0 [Ai(61 + 7ao - 20)]} . 



In the above expressions, P — fi^ — iJ,9 — hq (l — j i). 



(B23) 



(B24) 



(B25) 
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